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Abstract 

The structural behaviour of CsCdFs under pressure is investigated by means of theory and 
experiment. High-pressure powder x-ray diffraction experiments were performed up to a maximum 
pressure of 60 GPa using synchrotron radiation. The cubic Pm3m crystal symmetry persists 
throughout this pressure range. Theoretical calculations were carried out using the full-potential 
linear muffin-tin orbital method within the local density approximation and the generalized gradient 
approximation for exchange and correlation effects. The calculated ground state properties - the 
equilibrium lattice constant, bulk modulus and elastic constants - are in good agreement with 
experimental results. Under ambient conditions, CsCdFa is an indirect gap insulator with the gap 
increasing under pressure. 

PACS numbers: 62.50.-p, 61.05.cp, 62.20.de, 71.15.Nc 
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I. INTRODUCTION 



Ternary fluorides with the perovskite crystal structure have been extensively studied 
over several decades, as they have several potential applications because of their optical 
propertiesi*^, high-temperature super-ionic behaviour^, and physical properties, such as 
ferroelectricity^, antiferromagnetism^, and semiconductivity®. The application of CsCdFs in 
the field of luminesence^i^i^ has motivated several experimental and theoretical investigations 
of defect structures involving 3d transition metal ions.— 'ii'i^'^^i'i^ii^'i^ii^'i^. Elastic properties, 
including thermal expansion and acoustical measurements of higher order elastic constants, 
have been reported^i^SiSiiS^. Perovskite fluorides may exhibit structural phase transitions 
as function of temperature and pressure.— The present work is a combined theoretical 
and experimental study of the ground state and high-pressure properties of CsCdFa. We 
present the equation of state resulting from high-pressure diamond-anvil cell experiments 
up to 60 GPa. The ideal cubic structure is preserved over this entire pressure range. We 
also present the equation of state, the elastic constants and the electronic structure as ob- 
tained from theoretical calculations based on density functional theory using two different 
approximations for the exchange-correlation functional. 

The remainder of the paper is organized as follows. Details of the computational method 
as well as details of the experimental setup are outlined in section II. The measured and 
calculated equations of state are presented in section III together with calculated ground 
state properties and elastic properties. The electronic structure and the pressure variation 
of the band gap are discussed in section IV. Finally, conclusions are given in section V. 



II. COMPUTATIONAL AND EXPERIMENTAL DETAILS 



A. The electronic structure method 



The total energies and basic ground state properties of CsCdFa were calculated by the 
all-electron linear muffin-tin orbital (LMTO) method^^ in the full-potential implementation 



of Ref. l27|. In this method, the crystal volume is split into two regions: non-overlapping 
muffin-tin spheres surrounding each atom and the interstitial region between these spheres. 
We used a double n spdf LMTO basis to describe the valence bands, i. e. atom-centered 
Hankel functions with characteristic decay rate denoted by n are matched to a linear combi- 



nation of products of numerical radial function and spherical harmonics within the muffin-tin 
spheres. The calculations included the 5s, 5p, 6s, and 5d partial waves for cesium, the 5s, 
5p, and 4d partial waves for cadmium, and the 2s and 2p partial waves for fluorine. The ex- 
change correlation potential was calculated within the local density approximation (LDA)^ 
as well as the generalized gradient approximation (GGA) scheme^^. The charge density and 
potential inside the muffin-tin spheres were expanded in terms of spherical harmonics up 
to lmax=Q, while in the interstitial region, they were expanded in plane waves, with 14146 
waves (energy up to 156.30 Ry) included in the calculation. Total energies were calculated 
as a function of volume for a (16 16 16) k-mesh containing 165 k-points in the irreducible 
wedge of the Brillouin zone and were fitted to a second order Birch-Murnaghan equation of 
staters to obtain the ground state properties. 



The elastic constants were obtained 
conserving strains, as outlined in Refs. 



rom the variation of the total energy under volume- 
and 
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B. Experimental details 

Polycrystalline samples of CsCdFs were synthesized by the solid state reaction method as 
described elsewhere .-i^i^ For high pressure powder diffraction experiments, samples with a 
few ruby chips were loaded in a Mao-Bell type diamond anvil cell in a rhenium gasket (135 
/im hole diameter and pre-indentation thickness 65 /im). Pressure was generated with 325 
/im culet diamonds. Silicone fluid was the pressure transmitting medium.— X-ray diffraction 
experiments were performed at Sector 16 ID-B, HPCAT, at the Advanced Photon Source. 
The pressure inside the diamond anvil cell was determined by the standard ruby lumines- 
cence method.— The incident monochromatic x-ray wavelength was A = 0.36798 A. The 
x-ray diffraction patterns were recorded on an imaging plate with a typical exposure time 
of 10-20 s, with an incident beam size of 20 x 20 /im. The distance between the sample and 
the detector was calibrated using a Ce02 standard. The patterns were integrated using the 
Fit2D software program^'' and the cell parameters were obtained using the JADE package. 

The evolution of diffraction patterns as a function of pressure is shown in Fig. 1. At ambi- 
ent conditions the cell parameter a = 4.4669(7) A obtained experimentally agrees well with 
earlier reports provided in the Inorganic Crystal Structure Data base (ICSD).— Thermal 
expansion and temperature dependent structural experiments on CsCdFa were performed 



earlier by Reddy et al.— They found that the coefficient of thermal expansion increases 
with temperature. Rousseau et al.— investigated the crystal structure and elastic constants 
of RbCdFa, TlCdFa and CsCdFs using x-ray diffraction and Raman techniques at varying 
temperatures. They found tetragonal distortions in RbCdFs and TlCdFs at 124 K and 191 
K respectively. No structural phase transition were reported in CsCdFa in their temperature 
dependent study. Investigations by Haussuhl^ showed that temperature induced transitions 
can be observed in the pseudocubic perovskite CsCd(N02)3, which undergoes a phase tran- 
sition from R3 to PmSm symmetry around 464 K.^^ These reports indicate that the cubic 
structure of CsCdFa is relatively stable in comparison with RbCdFs, TlCdFa and KMgFa 
for a wide range of temperature from 4 K to 300 K. 

Fischer— investigated the pressure dependence of the second order elastic constants of 
CsCdFs and related fluoroperovskites. He reported a decrease in the force constant in 
RbCdFs, which is connected with the instability of the F ions. Such instability has not been 
noticed for CsCdFs, however a possible phase transition at high pressure was suggested.— So 
far, no high pressure experiments have been reported to ellucidate the structural behaviour 
of CsCdFa. In a previous work we showed that cubic KMgFa is stable up to 40 GPa^'^. If 
any structural phase transition is to be induced in CsCdFs, one would expect that pressures 
above 40 GPa are required. However, the diffraction patterns collected in our experiments 
show no pressure induced structural transformations up to 60 GPa. That is, similarly to 
KMgFs, the cubic phase of CsCdFa remains stable up to the highest pressures achieved in 
our experiments. 

The bulk modulus was obtained by fitting the pressure-volume data to a second order 
Birch-Murnaghan equation of state:— 

P = ^Bo [{V/Vo)-'/' - (V/Vo)-'/'] ■ {1 + 3/4(5^ - 4)[{V/Vo)-'/' - 1]}, (1) 

where Bq is the bulk modulus and B'q its pressure derivative. Least square fitting resulted 
in a Bq of 79(3) GPa, with a B'^ of 3.8. Thus, CsCdFs has a slightly higher bulk modulus 
than KMgFg ( Bo = 71.2 GPa, B'^ = 4.7)^ 
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FIG. 1: (Color online) Powder x-ray diffraction patterns of CsCdFa recorded at ten pressures up 
to 60 GPa 

III. GROUND STATE AND ELASTIC PROPERTIES 

Powder x-ray diffraction patterns collected at several pressures are shown in Figure 1. On 
compression, the diffraction patterns remain unchanged up to 60 GPa, except for the shifts 
of diffraction lines caused by the decreasing lattice constant and the broadening due to the 
intensity decrease with increasing pressure. This implies that no structural transformations 
occur up to 60 GPa in CsCdFs. Figure 2 shows the measured equation of state of CsCdFs 
and compares it with theoretical curves calculated within the LDA and GGA. As is typical, 
the LDA leads to slightly smaller volume at a given pressure than the experiment, while the 
GGA gives a larger volume. 

The lattice constant and bulk modulus measured in the present work as well as values 
calculated within the LDA and GGA are given in Table L Results from earlier experimental 
works are also quoted for comparison. The bulk modulus obtained in our experiments Bq = 
79(3) GPa with B'q = 3.8, compares well with values calculated using the LDA. The lattice 
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FIG. 2: (Color online) Measured (circles) and calculated (full line: LDA, broken line: GGA) 
pF-relation of CsCdFa. 

constant obtained within the LDA is L6 % lower than the experimental value, while the 
corresponding bulk modulus is 4.7% lower than the experimental value, which is the usual 
level of accuracy of LDA. When comparing the results obtained within GGA, the lattice 
constant is 2.2 % higher than the experimental value and the corresponding bulk modulus is 
32% lower than the experimental value. Since the present theory does not describe thermal 
effects, it is more relevant to compare to the low temperature experimental lattice constant, 
which is 0.4 % lower than the room temperature (RT) value (Table I). Further, correcting the 
low temperature experimental lattice constant for the zero-point motion of the ions, which 
is not considered in the theory, reduces the best value to compare to theory even more. 
Hence, the LDA lattice constant is in fact closer to experiment than the GGA value. The 
excellent agreement between the LDA calculated bulk modulus and the experimental value 
is, however, fortuitous, since the bulk modulus depends strongly on volume. Due to the 
overestimated equilibrium volume within GGA (and underestimated with LDA), an error 
is introduced in the calculated bulk modulus. Therefore, we recalculated the bulk modulus 
also at the experimental volume, using the simple scaling relation:— 




fi 



The corrected values are also quoted in Table 1. We find that this diminishes the discrepancies 
between the LDA and GGA results, as expected. In addition, the LDA bulk modulus now 
becomes sma//er than the GGA one for CsCdFs, and both functionals are seen to actually 
underestimate the bulk modulus, by approximately 25 % (LDA) and 7 % (GGA). 

The elastic constants of CsCdFa calculated within LDA and GGA are listed in Table 2 
where they are also compared to experimental results as well as earlier calculations. The 
LDA overestimates the Cn value by 28% and the corresponding GGA value is 1.8% lower 
than the experimental value. The C12 value within LDA is 5.4% lower than the experiment, 
and the corresponding GGA value is 33.3% higher than the experiment. The C44 elastic 
constants are higher by 9% both within LDA and GGA compared to experiment.— The 
Cii elastic constant obtained within GGA is much closer to the experimental value than 
the LDA value, while for Cu the situation is reversed, and for Cu the two approaches 
give approximately the same value, and in reasonable agreement with the experimental 
value. The elastic constants depend sensitively on the volume as illustrated in Figure 3 
for the case of LDA. Hence the above values carry an appreciable uncertainty inherited 
from the volume inaccuracy of the LDA or GGA approach. The pressure derivatives of the 
elastic constants are more stable quantities. They are compared to experimental values^^ 
in Table III with a reasonable agreement. The Cu has the strongest pressure dependence 
in accordance with the experiment, while the theoretical C44 value is almost constant with 
pressure, while experiment find it increasing with pressure. A point of caution is the fact 
that the present theory pertains to T = K, while experiments are performed at room 
temperature. Finite temperature generally tends to reduce the elastic constants because 
of thermal expansion. Using the calculated elastic constants we calculated the anisotropy 
factor A = 2Cu/iCu - C12). We find an A = 0.49 for LDA and A = 0.70 for GGA. The 
experimental^^ value calculated from the elastic constants is 0.74 which is slightly higher 
than the GGA value. 

IV. ELECTRONIC STRUCTURE 

The calculated electron band structure of CsCdFs is shown in Figure 4, while the density 
of states is displayed in Figure 5. The valence bands consist of the fluorine p bands, which 
dominate the first 2 eV below the valence band maximum. Below these bands appear the 
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FIG. 3: Volume variation of the clastic constants of CsCdFa as calculated in LDA. The volume 
is given relatively to the experimental volume, while the theoretical LDA equilibrium volume is 
marked with a vertical bar on the first axis. 



Cs 5p bands and the Cd 4d bands, which hybridize heavily in the region from 6 to 4 eV 
below the top of the valence band. The Cd 6s states dominate the bottom of the conduction 
band. The insulating gap is calculated to be 3.16 eV in LDA (3.67 eV in GGA). With 
compression, the gap increases almost linearly, at the rate 

The conduction band minimum occurs at the F point, while the valence band maximum 
occurs at the R-point (1/2,1/2,1/2)^, however almost no dispersion is found for the up- 
permost valence band along the line connecting the R-point with the M-point (1/2, 1/2, 0)^. 
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FIG. 4: Band structure of CsCdFa (using LDA, at the experimental lattice constant, a = 4.4669 
A). The zero of energy is set at the position of the valence band maximum. 

V. CONCLUSIONS 

In the present work, a combined theoretical and experimental analysis of the structural 
stability and equation of state of the fluoroperovskite CsCdFs has been carried out up to 
a pressure of 60 GPa. We find that the compound remains in the cubic structure in the 
entire pressure range studied. The calculated equilibrium lattice constant, bulk modulus 
and elastic constants agree well with available experimental data. Our electronic structure 
calculations show that CsCdFa is an insulator with an indirect gap which increases with 
pressure. 
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FIG. 5: (Color online) Valence band density of states of CsCdFs (using LDA, at the experimental 
lattice constant). The zero of energy is set at the position of the valence band maximum. The 
partial projections onto the spheres of Cs, Cd and F are shown with dotted (blue), dashed (red) 
and dash-dotted (green) lines, respectively, while the full line gives the total density of states. Units 
are electrons per eV and per formula unit. 
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TABLE I: Lattice constants (in A), bulk modulus Bq (in GPa) and its pressure derivative Bq, of 
CsCdFs as obtained by experiment and theory. The bulk moduli have been calculated both at the 
experimental and theoretical volumes {Bq{Vq^^) and Bq{Vq^^), respectively) 

Lattice constant ^o(V'o'') Bo{Vq'''^) B'^ 

GGA" 4.567 53.3 73.6 4.9 

LDA" 4.397 75.6 59.7 4.9 

Expt.(RT) 4.4669(7)'', 4.465^ 4.4662^ 79(3)" 3.8^,5.8°' 



Expt.(r = K) 4.452^ 
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TABLE IL Calculated elastic constants and shear modulus (G), all expressed in GPa, of CsCdFs 
at the theoretical equilibrium volume. 

Cii Ci2 C44 G 

GGA 105.8 27.0 27.7 32.4 Present 

LDA 150.2 38.3 27.5 38.9 Present 

Expt. 107.8±0.2 40.5±0.5 25.0±0.2 28.5±0.5 Ref. 19 



TABLE IIL Calculated LDA pressure derivatives of the elastic constants of CsCdFs at the exper- 
imental equilibrium volume. 
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